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ABSTRACT 



It has been noted for a long time that the spectra of observed continuum 
emissions in many solar flares are consistent with double power laws with a 
hardening at energies ~> 300 keV. It is now largely believed that at least in 
electron-dominated events the hardening in photon spectrum reflects an intrinsic 
hardening in the source electron spectrum. In this paper, we point out that a 
power law spectrum of electron with a hardening at high energies can be explained 
by diffusive shock acceleration of electrons at a termination shock with a finite 
width. Our suggestion is based on an early analytical work by Drury et al., where 
the steady state transport equation at a shock with a tanh profile was solved for 
a ^-independent diffusion coefficient. Numerical simulations with a p-dependent 
diffusion coefficient show hardenings in the accelerated electron spectrum which 
are comparable with observations. One necessary condition for our proposed 
scenario to work is that high energy electrons resonate with the inertial range 
of the MHD turbulence and low energy electrons resonate with the dissipation 
range of the MHD turbulence at the acceleration site, and the spectrum of the 
dissipation range ~ k~^''^. A ~ /c"^'^ dissipation range spectrum is consistent 
with recent solar wind observations. 

Subject headings: Sun: flares — Sun: particle emission — Sun: X-rays, gamma rays 
— acceleration of particles 
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Introduction 



14 Our Sun is an efficient particle accelerator. Ions with energy up to ~ Gel^/nucleon are 

15 detected in-situ during large Solar Energetic Particle (SEP) events where both large flares 

16 and fast CMEs often occur together. At flares, el ectron Bremsstrahlun g is believed to b e 

17 the main source of the continuum radiation (e. g. 



Ramaty et al. 



(119751 ): IVestrandl (119881 ) I 



18 Continuum emissions provide invaluable information that constrains the underlying 

19 acceleration mechanism. These constraints include the energy budget, the t otal n umber 



20 of electrons, th e acce leration time scales, etc. See reviews by 



Miller et al 



Zharkova et al. 



(Il997h and 



(120111 ) for a detail discussion on various acceleration mechanisms and the 



22 implications of these constaints on them. 

23 One observational constraint that received less attention, even though has been noted 

24 for a long time , is the hardening of the continuum spectrum at high energies (often ~> 300 



keV). In 1975, 



Suri et al 



(Il975[ l examined the X-ray and gamma-ray flux in the August 
4, 1972 event and concluded that the X and gamma-ray flux was produced by a single 
population of electrons with a break in its sp e ctrum , instead of two separate populations 



acting independently. Later, 



Yoshimori et al. 



( 119851 ) , using Hinotori spacecraft, examined 



the hard X-ray (HXR) spectrum in a broad energy range (20 keV - 7 MeV) for four flares 
that showed significant h arden ing at energies above ~ 400 keV. They confirmed the earlier 



31 suggestion of 



Suri et al 



( Il975l ) that the hardening in the continuum reflects an underlying 
hardening in the source electron spectrum. Spectral hardening also occur in events where 
there are clear signatures of gar nma-ray lines. The most recent report of such an event is 



from the Fermi observation (see 



Ackermann et al. 



(120121 )) where spectral hardening was 



35 found at above several hundred keV. 

36 Note that hardening in the source electron spectrum is not the only cause for 

37 a hardening in the photon spectrum. Various processes, such as electron-electron 
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BremsstrahlunK , positronium annihilation continuum and inverse 



38 Bremsstrahlung, proton 

39 Compton emissions (see IVestrandl (jl988l )). may all lead to some hardening of a continuum 

40 emission from a straight power law (without hardening) source electron spectrum. However, 

41 for parameters appropriate to a solar flare site, the contributions of these processes 



are relatively smal 



43 ( Kontar et al. 



and the resulting hardening of the spectral index is perhaps ~ 0.5 



20071 ). Therefore the se processes can n ot explain events where the change 



44 of spectral indices are as high as 2. 



Park et al 



(119971 ) studied photon spectral hardening 



around 1 MeV for four flares. In their scenario, the emission is a simple sum of the thin 
target emission from the trapped electrons at the acceleration site near the loop top and the 
thick target emission from the escaping electrons precipitating on the solar surface. With 
the assumption that electrons having smaller energies have shorter escape times, and noted 
that the energy dependence of the Bremsstrah lung cross section differs in the nonrelativistic 



and the relativistic regimes. 



Park et al 



51 hardening. Note in the scenario of 



19971) were able to account for the observed spectral 



Park et al. 



( Il997l ). the energy dependence of the escaping 



52 time decides the hardening, and the accelerated electron spectrum in t he accelerated regi on 



Moses et al. 



f ll989h 



does not (need to) have a spectral hardening. In-situ observations by 
however, showed that electron spectral hardening is rather common in short duration events. 
Assuming these in-situ electrons are the source electrons es caping from the acceleration site 



56 through interchange reconnection, then 



Moses et al 



(jl989l )'s results also suggest that the 



accelerated electron population has a hardening at high energies. 

High energy electrons also lead to microwave e r nissio ns through gyro-synchrotron 



59 radiation (see e.g. the recent review by 



White et al. 



(120111 ) for a detailed discussion of the 
relationship between solar radio and HXR emissions). In an early study , usin g BATSE 



Silva et al. 



( I2OOOI ) examined 



61 (HXRs) and Owens Valley Radio Observatory (microwaves) 

62 27 solar flares with multiple peaks (a total of 57) which were observed at both HXR and 

63 microwave wavelengths. Fitting the HXR spectra by a single power law and the microwave 
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64 spectra as gyrosynchrotron emissions, ISilva et al.l (|2000( ) found that in 75% of the bursts, 

65 the inferred spectral indices of the electron energy distribution of the microwave-emitting 

66 electrons w e re ha rder (by 0.5 to 2.0) than those of the lower energy HXR emitting electrons. 



Silva et al. 



( I2OOOI ) concluded that there exists a breakup in the energy spectra of the source 

68 electrons at around ~ 300 keV, in agreement with previous observations of HXR-alone 

69 spectra of giant flares. 

70 Note, however, in most events the HXRs are emanated from the footpoints of flare 

71 loops and microwaves are emanated from the tops of flare loops. Furthermore, there is also 

72 a delay between the peak of the microwave emission and the HXR emission. So there are 

73 transp ort effects between electrons generating microwave emissions and those generating 



74 HXRs 



White et al 



( 120111 ). Both the harder HXR spectrum at the footpoints and the delays 

75 of the microwave emission could be caused by magnetic trapping of higher energy electrons 

76 near the loop top and the precipitation of lower energy electro n to the footpoints, a s first 



77 suggested by 



Melrose and Brown 



Kawate et al 



(12012), HXR 



(119761 ). In a recent study by 

78 and microwave emissions from 10 flares were analyzed. Although the emissions were at 

79 different locations and the spectral indices for microwave emissions are harder than those of 

erated electrons that is consistent with the 



Kawate et al. 



( I2OI2I ) were able to produce 



80 the HXRs, by assuming a spectrum for the acce 

81 HXR emissions (but extend to higher energies), 

82 microwave spectra comparable to the observations. The authors concluded that it is a 

83 single electron population that is responsible for the HXRs and microwaves emissions and 

84 the hardening of the microwave emission is due to a more efficient trapping of electrons 

85 with higher energies. In another study, to minimize the effect of the tr apping of h i gh 



Asai et al. 



(120 131 ) 



86 energy electrons on the resulting spectra of looptop microwave emissions 

87 examined both the HXR and microwave spectra prior to the peak emission in 12 flares. 

88 They still find a significant hardening of the source electrons for the microwave emissions. 

89 These authors suggest that there is an intrinsic spectral hardening for the source electron 
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90 spectrum around several hundreds of keV and the microwave gyrosynchrotron emission is 

91 due to electrons at higher energies (in the harder part of the spectrum). 

In this wo r k, we do not consider microwave emissions and focus on HXRs alone. 



Vestrand et al. 



(119991 ) identified 258 flare events using SMM observations. Among these, 



94 many are electron-dominated events with no clear sign ature of gamma-ray emissions 



95 (IRieger fc Marschhauser 



1990 



Marschhauser et al 



1994 ). In these events, the contribution 



96 of nuclear gamma ray lines is minimal and the continuum is mainly due to Bremsstrahlung 

97 of the energetic electrons. The spectral hardening can be clearly seen in many of those 

98 electron-dominate events. A careful examination of these events based on the mechanism 

99 proposed here will be reported elsewhere (Kong et al. to be submitted). 

100 A hardening in the source electron spectra is hard to explain for any acceleration 

101 mechanism. In this paper, we propose a scenario which is based on diffusive shock 

102 acceleration (DSA) to explain the observed spectral hardening. 



Electron acce 



eratio n at a termination shock (TS) in solar flares is not a new idea. 



Tsuneta &: Naitd (11998) wer e the first to consider electron acceleration via DSA at a flare 



105 TS. iTsuneta fc Naitol (jl998l ) pointed out that slow shocks bounding the reconnection 

106 X-point can heat the plasma up to perhaps 10-20 MK, providing abundant seed population 

107 which are accelerated to 1 MeV in 0.3 to 0.6 seconds at the TS. 



Mann et al 



(12009 ) 



Mann et al. 



108 Noting that the standing TS is of quasi-perpendicular in nature, 

109 consid ered shock drift acceleration (SDA) at a standing TS. In the work of 
no (120091 ). most energy gain of electrons is through a single reflection at the shock front. 

111 Therefore to accelerate electrons to high energies, a stringent requirement of 9bn (e.g. 

112 6bn > 88°) is needed. However, while the TS on large scale is of quasi-perpendicular, 

113 small scale structures (such as ripples) exist on the shock front. Indeed, the plasma in the 

114 reconnection region is unlikely to be homogeneous, so the resulting TS is unlikely to be 
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planar. Recently, iGuo fc Giacalond ( 120121 ) ha ve examined electron acc eleration at a flare 



Guo fc Giacalond ( 120121 ). many small scale 



TS using a hybrid code. In the simulation of 
ripples were identifled along the shock surface. The existence of these ripples suggests 
that assuming a shock with a. 6bn > 88° across the shock surface may be unrealistic. 
Furthermore, the existence of these small scale structures implies that one single fleld line 
can intersect the shock surface multiple times. Co nsequently, the accelerat ion process will 



Tsuneta fc Naitd (Il998[ ) and assume 



be of diffusive in nature. In this work, we follow 
the electron acceleration at a flare shock can be described by the DSA mechanism. Note 
that the existence of a TS in a flare site is not trivial. Observational evidence of flare 
TS has been reported by 



Warmuth et al. 



( l2009l ). who used dynamic radio spectrum from 



the Tremsdorf radiospectrograph to show that there was a type-II radio bursts from a 
standing TS at ~ 300 MHz during the impulsive phase of the XI. 7 flare of 2001 March 
29. Besides the existence of a TS, the area of the TS shock has also to be large enough 
(~ 10^° cm^ in large flares ) to account for the observed flux of HXRs generated by high 
ene rgy electrons. By ass uming a 50% contour of the N RH source at 327 MHz (see Figure 2 



130 of ( IWarmuth et al. 



20091 1 ) being a proxy for the shock. 



Warmuth et al. 



( I2OO9I ) estimated a 



shock area of ~ 1.3*10^° cm^ in the 2001 March 29 XI. 7 flare. We do note that these areas 
are much larger than the areas of HXR sources and are more comparable to active region 
sizes. Whether or n ot the size of the TS ca n be this large remains to be examined. Using 



134 the same technique. 



Warmuth et al. 



( 120091 ) nevertheless obtained similar areas for other 
events where TS were observed. Of course, for smaller flares (like M flares), the area of the 
active region is smaller and we expect the area of the shock is also smaller. 

Besides shock acceleration, models based on stocha stic (aka 2nd-order Fermi) 



138 acceleration exist. For example. 



Miller et al. 



(11996 



19971 ) assumed the presence of some 



139 large scale turbulence at the flare s ite and considered the co upled system of the wave 



cascading and particle acceleration. 



Miller et al 



(11996 



19971 ) showed that various modes 
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141 of waves (Alfvenic and fast mode waves), as they cascade to small scales, can efficiently 



142 accelerate both i ons a n d electrons. Simila r proces ses have also been studied 



Petrosian et al 



(19941 



Park et al. 



Park et al. 



(119971 ). Unlike 



Miller et al 



( 11996 



19971 ) 



jy, for example, 



Petrosian et al. 



(119971 ) did not address the cascading of the turbulence and assumed the 



145 wave spectra is given. 



In this work, we do not con si der stochastic acc eleration. However, as in 



( 11994h 



Miller et al 



( 11996 



19971 ): 



Park et al. 



Petrosian et al. 



(Il997l ). we assume the diffusion coefficient k is 



decided by the underlying turbulence power at a flare site. 



Diffusive shock acceleration of electrons at a finite-width termination shock 



At a piecewise shock, the standard steady state DSA predicts a power law spectrum 
~ for energetic particles. The power law spectral index a is given by 3s/(s — 1), where 
s = Ui/u2 is the compression ratio, Ui and U2 the upstream and down stream flow speed i n 



Drury et al. 



(119821 ) 



the shock frame. In the case of a shock having a flnite width ~ L^iff, 
showed that the spectral index depends on the shock width. Assuming the background fluid 
speed is given by a tanh proflle: 

Ui + U2 Ui — U2 



U[X] 



-tanh{x / 



diff) 



1 



then the spectral index a becomes (jPrury et al 



1982|), 



3s 



■(1 



1 1 



where (3s is a dimensionless parameter and is related to the diffusion coefficient k through 



(2) 



K = (3s{u - Ui){u - U2) 



dx 
du 



Ml - U2 



L 



diff- 



(3) 



150 Although 



Drury et al. 



( 119821 ) considered only the case of p-independent k where analytical 



151 solutions can be obtained, one can see from the above that for a k increasing with p the 



- 9 - 



152 spectrum will harden at high energies. Because the factor of in equation ([2]), the 

153 spectral index quickly approaches the limit of 3s/(s — 1) when (5s > 1. When [5^ is small, 

154 however, the second term in the bracket of equation ([2]) dominates and the spectrum can 

155 be very soft. 

156 Clearly the momentum dependence of the diffusion coefficient n decides the shape of 

157 the spectrum. At a flare site, the /t of energetic electrons is decided by the turbulence 
15B level. At large scales, the turbulence is of Alfvenic and particle-wave cyclotron resonance 
159 can ac celerate ions to high energies via the stochastic acceleration process (e.g. 



Miller et al, 



160 f ll997l )). For electrons, except at very high energies, however, they do not resonate with 



161 Alfven waves, t 



lerefore they inte ract with other waves, for example, fast mode and/or 



162 whistler waves ( iMiller et al. 



19961 ). 



Note that 



Drury et al. 



(Il982[ ) did not consider the effect of the energetic electrons 

164 on the shock. In a more refined and self-consistent analysis, the pressure of the energetic 

165 electrons needs to be taken into account and it will affect the shock width. This is 

166 similar to a modified shock struct ure caused by energetic cosmic rays as first examined by 



Axford. Leer fc McKenzie 



(Il982l ). Such a discussion, however, exceeds the scope of this 



168 work and we do not consider the back reaction of energetic electrons on the shock structure. 

We assume the turbulence at a flare, as in the solar wind, is described by an inertial 
range joining to a dissipation range and the power density I{k) is given by. 



i{k) = i{ko){{lr^H{h -k) + {^r^H{k - k,)), 



(4) 



where and are the spectral indices in the inertial range and dissipation range, 
respectively. We assume = 2.7 (see below) and consider three cases for ej: 5/3, 1.5 and 
1.0. The case of = 5/3 corresponds to a Kolmogorov cascading; the case of = 1.5 
corresponds to a Iroshnikov-Kraichnan (IK) cascading, and the case of = 1.0 corresponds 
to a Bohm-like diffusion (see below). At very small k, the energy containing range sets in 
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and I{k) bent over. The normalization of I{k) is given by, 

/ + 00 
I{k) =< 6B^ > . (5) 
-oo 

For a wide range of electron energy, the resonating wavenumber k is in the dissipation 
range. In the solar wind, one finds a spectrum ^ /c~^'^ to ^ ° i i i the dissipation 



range (ILeamon et al 



1998 



1999 



Chen et al. 



2010 



Howes et al. 



2011 



Alexandrova et al 



2OI2I ). Unlike the inertial range, the nature of the turbulence in the dissipation range 



is still unde r debate. Two possible sc e narios include Landau dam ping of kinetic Alfven 



waves e.g. 



Leamon et al. 



( IStawicki et al. 



2001 



1999 



2000 



Bo 



Krishan fc Mahajan 



dvrey and Perez 



2004; 



Galtier 



20121) , or whistler waves e.g. 



20061 ) . For KAWs, k^ > k\\ 



electron-wave interaction is through the Landau resonance and KAWs can effectively heat 
electrons. Whistler waves have Up < u < Q and electrons can interact with whistler waves 
through the cyclotron resonance. The resonance condition is. 



u — k\\V\\ 



(6) 



where VL = eB/i^nie) is the electron cyclotron frequency and 7 is the Lorentz factor. For 
low frequency waves u <VL, the resonance condition (on taking n = 1) yields 



= ^/\k\\ 



(7) 



where /i is the pitch angle of the electron. Note from equation ([6]), one can see that when 
the energy of electron is high enough, it can also resonate with Alfven waves. In this work 
we assume the dissipation range turbulence is whistler-wav e-like and elec t rons can resonate 



with the wave t 



(I2OO3I) 



i rough the cyclotron resonance. As done in 



Li et al. 



Gordon et al. 



fll999[ ): 



Rice et al. 



( I2OO5I ) , we further simplify the resonance condition by replacing k = Vt/ jiv 
with k = Q/v, which corresponds to an extreme resonance br oadening. Th e pitch angle 
diffusion coefficient D^^, from the Quasi-linear Theory (QLT) (jjokipiil 119661 ) is, 

1 - /i2 fi2 



D 



MM 



\n\v Bl 



i{k = n/i2v) 
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The diffusion coefficent k is related to -D^^ through, 

•2 '•+I(l-/i2)2 



V 



K 



1 



D 



lQ9?I{k = n/v) 



(9) 



We make no attempts to estimate the turbulence level at the reconnection site in this 
work. Instead, we are more interested in the energy dependence of k. From equation (Q, 
we have 



K = Kq- 



{p/Po 



7 



(10) 



169 where subscripts i or d denote whether electrons resonate with the inertial or the dissipation 

170 range of the turbulence. For electrons resonating with the dissipation range that have 

171 a ed ~ 2.7, equation ( ITOl) suggests that k has a very shallow dependence on electron 



172 momentum (energy). In comparison, for electrons reso nating with the inertia 



Tsuneta fc Naitol f ll998f ). the Bohm 



range, k 



173 increases quickly with particle momentum (energy). In 

174 diffusion approximation was used, in which case k ~ vRi, where Ri is electron's gyroradius. 

175 This corresponds to an = 1. 

The fact that k has a very shallow dependence on the electron's momentum in the 
dissipation range and a strong dependence in the inertial range is the key to understand the 
hardening of electron spectrum. In Figured] we plot Ps as defined in equation ([2]), where 
from equation ffTOj) . we have 



/3s = /3o 



{p/Po 



7 



(11) 



176 We set Po = 0.2. This value yields an electron spectral index at low energy to be ~ p 

177 comparable to flare observations. 

Note, from eq u ation ([3]), Ps also depends on the width of the shock. Simulations by 



Scholer fc Burgess 



( 120061 ) suggested that the shock width is of the order of ion inertial 

180 scale length ~ (c/upi). On the other hand, observations of the Earth's bow shock (at 

181 quasi-perpendicular configurations) showed that its ramp width is somewhat smaller than 
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Newburv et al 



Zank et al 



c/cup,) (jScudder et al.lll986l : iBalikhin et al.lllQQSt iNewbury et al.lll998[ ). In particular 



■998[ ) found considerable fine structures of the order of ~ (c/wpe) 



( 120011 ) suggested that these fine structures will help to circumvent the injection 



185 problem for Anomalo us cosmic rays. In a very recent study, using Clusters observation, 



Schwartz et al. 



(120111 ) showed that at the Earth's bow shock half of the t emperature 



Schwartz et al. 



187 occurr ed in about ~ Ic/ujpe or ~ {l/l)c/ujpi. The total width of the shock in 

188 (120 111 ), which is close to Ldiff in our work, however, is another factor of ~ 6 (see their 

189 figure (3)). Therefore, in this work, we assume the shock width is given by the ion inertial 

190 length scale L^iff ~ c/cUpj. 

191 The break point ph in Figure [1] is ph ~ l^e^/kb with the wave number separating 

192 the inertial range and the dissipation range. The scale at which the inertial range transits 

193 into the dissipation range is still a much debated iss ue. It has been argued th at it could 



194 be the thermal proton Larmor radius 



(ILeamon et al 



195 


ion inertial 


ength 


196 


Smith et al. 


2001) 



ength ~ ^ with f2p the proton cy clotron frequency ( Leamon et al. 



200ll ). Consider a typical flare site (IMiller et al. 



1996 



1998 



1999h or the 



Mann et al. 



2000 



120091 ) with 

197 a temperature of T ~ 5 MK, a magnetic field of i? ~ 200 Gauss and a density of rze ~ 10^ 

198 cm^'^, we find an Alfven speed Va ~ 1-38 * 10^ km thermal proton speed Vth ~ 200 

199 km proton gyrofrequency fip = 1.91 * 10^ Hz. Consequently, the thermal ion 

200 Larmor radius is ~ 0.10 m and the ion inertial length is ~ 7.2 m. If k^, is the reciprocal of 

201 the thermal ion Larmor radius, then p}^ ~ 0.64 MeV/c and the corresponding kinetic energy 

202 is 0.31 MeV. This is in good agreement to the observed continuum emission break locations 

203 ~ 300 keV. On the other hand, if fcf, is the reciprocal of the ion inertial length, then p\^ ~ 39 

204 MeV/c, much too hi gh for the proposed s c enari o. Therefore our proposed scenario favors 



205 the the suggestion of 



Leamon et al. 



(11998 



19991 ) that the dissipation range sets in at the 



206 thermal ion Larmor radius scale. Comparing to the width of the shock, which is the ion 

207 inertial length scale 7.2 m, the gyro-radius of an electron Ri = 'jv/Qe is 0.1 (0.4) m for a 
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20B kinetic energy of 300 keV (2 MeV). 

209 Using a momentum dependent (5s as in equation fITT]) . we numerically solve the 

210 steady-state transport equation. We set po to be 32 keV/c, which corre s ponds to an 



Drury et al. 



fll982h . given 



211 injection energy of 1 keV. We use the same shock profile as 

212 by equation ([T]) and assume a compression ratio of 3.5 (thus a strong shock). Note that 

213 the outflow plasma speed at a reconnection site is ~ Va- Therefore for a shock with 

214 a compression ra. t io of 3.5, ui — U2 in equation ([3]) is ~ 10^ km/s. We use /3o = 0.2. 



Tsuneta fc Naitol (11998 



a) have used the Bohm approximation for k. With the Bohm 

216 approximation k = ^vRi and the above values for a typical flare site, then a 300 keV 

217 electron will have /S^ = 0.2 and a 2 MeV electron will have (3s = 1.0, suggesting our choice 

218 of (3o = 0.2 is reasonable. 

219 Figure |2] plots the steady state electron spectrum for three cases that have different 

220 inertial range turbulence spectrum: i): Kolmogorov-like; ii): IK-like; and iii): Bohm 

221 diffusion approximation. In each panel, the two dashed lines are power law fittings 

222 ~ (p/po)~") with «! and a2 the fitted spectral indices to the spectrum at the low and high 

223 energies respectively, and Pm the fitted break momentum. We set pb to be ISpo ~ 0.416 

224 MeV/c in the simulation, which corresponds to a Eb of 0.15 MeV. Note that pm is larger 

225 than Pb = '-frrieQ/kb by about a factor of ~ 2. Figure |2] is the most important result of this 

226 paper. It shows that diffusive shock acceleration at a finite-width termination shock in solar 

227 flares can naturally lead to a hardening of the accelerated electron spectrum. 



22B 3. Discussions and conclusions 

229 Clearly, the hardening requires the following conditions to be met. First, the existence 

230 of a termination shock at flare site with a finite shock width L^iff ~ c/upi. Second, the 
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231 diffusion coefficient k needs to be close to a constant at low energies and increases with 

232 electron energy at high energies. Third, it is necessary that k < AULdiff at energies below 

233 the break and k > AULfUff at energies above the break. 

234 For any given ffare, none of these conditions are necessarily satisfied. 

235 Consider the first condition. While it is hard to identify a termination shock at a flare 

236 observationally, there are indirect clues of such shocks. For exam ple, type II radio bursts 

237 without frequency drift has been used by 



Warmuth et al. 



(120091 ) to infer the existence of 



238 flare termination shocks. Further observational evidence of flare termination shock, and in 

239 particular its size, are welcomed. 

240 For the second condition: if electrons resonate with the dissipation range of the 

241 turbulence through cyclotron resonance and that the dissipation range has a power 

242 spectrum I{k) ~ /c"^ '', then we find k, is indeed close to a constant at low energies and 

243 increases with electron energy at high energies. While in-situ solar wind observations do 

244 suggest such a ~ fc"^'^ dissipation range, direct confirmation of such a k dependence in the 

245 flare site is impossible. 

246 To satisfy the third condition will put a strong constraint on the turbulence level at the 

247 flare, which can vary much from one to another. Consequently, the hardening does not occur 

248 in all flares. For example, if for a given flare, /^^ ~ 1 instead of /3s ~ 0.2 at lower energies, 

249 then there will be no hardenings, even if both conditions 1 and 2 are satisfied. Because 

250 a larger Ps implies a larger k, therefore a less efficient acceleration, so one implication of 

251 our proposed scenario is the following: Events where the continuum emission extend to 

252 very high energies (efficient acceleration) likely show hardenings and have softer spectra at 

253 lower energies, and events where the continuum emissions do not extend to high energies 

254 (inefficient acceleration) likely have harder spectra at low energies than those events that 

255 extend to higher energies. 
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256 Another consequence of our proposal is the correlation between the low energy photon 

257 spectral index 71 and the break momentum p^. Consider two nearly identical flares A and 
25B B except that flare A has a larger ki, (i.e. the inertial range in flare A extends to a smaller 

259 scale). Then /3s at p < pm for flare A is smaller. Therefore, pm and ai are anti-correlated. 

260 Observations do show such a anti-correlation and this is discussed in details in Kong et al., 

261 (in preparation). 

262 In summary, we offer an explanation for the observed continuum spectral hardening 

263 in solar flares that is based on DSA. To our knowledge, no previous works have addressed 

264 the hardening of emission spectrum explicitly. Further observational and theoretical studies 

265 along the proposed mechanism will be pursued in future works. 
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267 grants NNH07ZDA001N-HGI and NNX11A064G at UAHuntsville and the 973 Program 
26B No. 2012CB825601 and NNSFC Grant Nos. 41274175 and 41028004 at SDUWH. XLK 
269 acknowledges financial support by the Shandong University Graduate Study Abroad Fund. 
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Fig. 1. — /3s as defined in equation ([3]). At low energies, electrons resonate with the 
dissipation range and Ps has a weak momentum dependence. Above pb, electrons resonate 
with the inertial range and Ps quickly increases with momentum. Three cases for the inertial 
range are considered. These are, from top to bottom, Bohm-like, IK-like and Kolmogorov- 
like. 



-17- 



lOkeV O.BBMeV lOMeV 1 OkeV 0.65WeV lOMeV 1 OkeV 55UeV lOMeV 



10= 


\ 




1 




10= 


t\ 








10= 




1 






10° 










10° 










10° 










10"= 
10-'° 


s 


»v \ 
S \ 

*^ V 

*s \ 






10"= 
10-'° 


i. 


V \ 

*v \ 
% \ 






10-= 
10-'° 




s \ 
^ \ 
X \ 






10-' = 










10-' = 






\ 




10-'= 






\ 

\ 




10"™ 


a, 


= 10.0 
= 5.6 


\ 

\ 




10"=° 


a, 


= 10,0 
= 5.2 


\ ^ 
\ 

\ 




10-™ 


a, 


= 10,0 
= 4.6 


\ 

\ 

\ 






Pm 


= 33.5pQ 


% 






Pm 


= 32.5p(, 


\ 






Pm 


= 29.5p(, 


\ 




10-=^= 

1 










,\ , ,„ 


10-== 








,\ , ,„ 


10-== 














10 


100 


1000 




10 


100 


1000 




10 


100 


10 



P/Pft P/Po P/Po 



Fig. 2. — The electron spectrum for three cases: i: Kolmogorov-hke inertial range; ii: IK-hke 
inertial range; iii: Bohm diffusion Approximation. 
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